Kondo effect in a one-electron double quantum dot: Oscillations of the Kondo current 

in a weak magnetic field 
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We present transport measurements of the Kondo effect in a double quantum dot charged with 
only one or two electrons, respectively. For the one electron case we observe a surprising quasi- 
periodic oscillation of the Kondo conductance as a function of a small perpendicular magnetic field 
\B±\ ;$ 50 mT. We discuss possible explanations of this effect and interpret it by means of a fine 
tuning of the energy mismatch of the single dot levels of the two quantum dots. The observed degree 
of control implies important consequences for applications in quantum information processing. 

PACS numbers: 72.15.Qm, 73. 21. La 73.23.Hk, 



The Kondo effect describes a bound state formed by 
interactions between a localized magnetic impurity and 
itinerant conduction band electrons shielding the local- 
ized spin. This results in an increased density of localized 
states at the Fermi energy, causing anomalous low tem- 
perature properties. In case of a degenerate ground state 
of a quantum dot (QD), the Kondo effect manifests itself 
as an enhanced conductance within the Coulomb block- 
ade region [H, d, This was first observed on large 
QDs with half integer spin 0, and later, for a total 
spin of S = 1, where the triplet states of a QD are de- 
generate [1, 0, HI- On a double quantum dot (DQD) a 
two- impurity Kondo effect was studied 

In this article we present the results of Kondo effect dif- 
ferential conductance (KDC) measurements on a DQD 
charged with one or two electrons in a perpendicular 
magnetic field B±. For only one electron (TV = I) in the 
DQD we observe a quasi-periodic structure of the KDC 
with a characteristic scale of B$ ~ 10 mT. In contrast for 
TV = 2 the KDC is found to be a monotonic function of 
Bj_. We discuss possible explanations for this effect that 
imply consequences in quantum information processing. 

Our sample is fabricated from an AlGaAs/GaAs he- 
terostructure. It embeds a two-dimensional electron sys- 
tem (2DES) with carrier density n s ~ 1.8 x 10 15 m -2 and 
electron mobility fi ~ 75m 2 /Vs (at T = 4.2 K) 120 nm 
below its surface. Figurc[ljb) shows Ti/Au-gates created 
by electron beam lithography. They are used to locally 
deplete the 2DES to define a one electron QD. The gate 
design is optimized for transport measurements through 
a QD charged by only few electrons [l(j. By decreasing 
the voltages applied to gates gc and gx (with respect to 
the 2DES) while increasing the voltages on the side gates 
gL and gR we deform the QD into a DQD (sketched in 
Fig. Efb)) [U The DQD is tuned to the regime 



of strong coupling to the leads and an order of magni- 
tude stronger interdot tunnel coupling of 2t n ~ 240 fxeV 
between the adjacent QDs [ll[. Measurements are per- 
formed in a dilution refrigerator at an electron tempera- 
ture T 2 des ~ 0.1 K. 

A nearby quantum point contact (QPC) is used to 
detect the charge distribution of the DQD shown in 
the stability diagram in Fig. QTa) It displays a 

lock-in measurement of the differential transconductance 
Gqpc — d/Qpc/dC/g L as a function of the dc voltages 
applied to gates gL and gR. In the lower left corner re- 
gion in FigJTJa) the DQD is uncharged (compare figure 
caption) [ll| . 

The differential conductance of the DQD is plotted in 
Fig. [He) as a function of the applied bias voltage Usu and 
the center gate voltage U gc . The DQD is tuned such, that 
the variation in U gc (x-axis) causes a shift in the stabil- 
ity diagram approximately along the arrow in Fig. QTa). 
Hence, a charge between TV = and 3 electrons, marked 
in Fig. []Jc) by numerals, is distributed symmetrically 
between the adjacent QDs. Within the diamond-shaped 
regions transport is impeded by Coulomb blockade (CB). 
Nevertheless, strong coupling of our DQD to its leads al- 
lows for inelastic co-tunneling causing an enhanced con- 
ductance within the CB regions at Usd ^ 0, (e.g. for 
TV = 1 at |e?7sD| > 2i ~ 240 ,ueV) 0. 

In addition, at Usu — an increased differential con- 
ductance is visible in the CB regions for TV = 1, 2 or 
3. We assign this zero bias anomaly to the Kondo effect 
on a DQD, here charged with only a few electrons. The 
observed KDC is small compared to the unitary limit 
(G -C 2e 2 /h). This is due to the tunnel barrier hindering 
electron transport between the two adjacent QDs and to 
an asymmetric coupling to the leads [HI]. For N = 1 
or 3 the KDC of the DQD can be described by the spin 
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FIG. 1: (Color online) (a) Stability diagram of the DQD. 
Plotted is the transconductance d/Qpc/d[/ gL (color scale) of 
the DQD as a function of the side gate voltages U gL and 
Ug n . The measurement was done by using a nearby QPC. 
A background is subtracted for clarity. Numerals denote the 
number of electrons charging the (left/right) QD. (b) SEM 
micrograph of the top gates on the sample surface. Arrows 
mark possible current paths through the DQD and the nearby 
QPC. The estimated DQD geometry is sketched in white, 
(c) Differential conductance dl/dUso of the symmetrically 
charged DQD as a function of bias voltage C/sd and center gate 
voltage U gc . The variation of U sc approximately corresponds 
to the path in the stability diagram indicated by an arrow in 
(a). Numerals indicate the number of electrons charging the 
DQD. 



1/2 Kondo effect, but for N = 2 the threefold degenerate 
triplet states lead to the KDC. This suggests that the 
exchange coupling separating the triplet states from the 
singlet ground state is smaller than either the Kondo or 
the electron temperature. 

Figured] displays the KDC at [/ S d ^ of the DQD as 
a function of a magnetic field B± perpendicular to the 
2DES for N = 1 and N = 2, respectively. Each point 
corresponds to the maximum KDC near zero bias mea- 
sured at constant gate voltage approximately along the 
white vertical lines in Fig. [ljc). Three of these traces 
G(Usd) are plotted in the inset of Fig. [2] For increas- 
ing B± the KDC is expected to monotonically decrease 
as the spin degeneracy is lifted. A theory by Pustil- 
nik and Glazman provides analytical expressions for the 
limits B <C -Bk and B Bk [ill , where the charac- 
teristic field Bk is determined by the Kondo tempcra- 
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FIG. 2: Kondo effect differential conductance in a one (two) 
electron DQD. The inset displays exemplary raw data curves 
of the KDC of the DQD at N — 2 as a function of the bias 
voltage for different magnetic fields but constant gate volt- 
ages. All raw data curves are measured within CB regions as 
sketched by the vertical white lines in Fig. [T] Black (N = 1) 
and grey (N = 2) circles in the main figure display the KDC 
in a logarithmic scale at Usd — (local maxima of raw data 
curves) as a function of a perpendicular magnetic field. Lines 
are model curves explained in the main text. 



For B± <C Bk, the KDC is de- 
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ture k B T K = g^ B B K - 

scribed by G ~ G (l - (B ± /B K ) 2 ) and for B ± > B K by 
G ~ Goo/ In 2 (Bx/B K ). G is the KDC at B ± = 0. The 
lines in Fig. [5] model these expressions with Tk = 0.1 K 
for N = 1 and T K = 0.12 K for N = 2. T K is taken iden- 
tical for both limits (solid and dashed lines), respectively. 
Being close to the electron temperature of the 2DES, Tk 
cannot be extracted from temperature dependences as 
usual. Nevertheless, the model curves used here are ex- 
pected to hold even for Tk ~ ?2des ■ The agreement with 
our data is satisfactory. 

For B± > 0.5 T the decrease of the KDC gets steeper 
due to a B± dependent decrease of the interdot tunnel 
coupling, specifically investigated for our DQD jll] . Tak- 
ing such effects into account does not change the Kondo 
temperature too much compared to the simple model pre- 
sented here. However, our simplified model causes the 
fit-parameter Goo to be strongly suppressed compared to 
G . 

Figure [3] shows detailed measurements of the KDC for 
N = 1 (main figure) and N = 2 (inset), as plotted in 
Fig. [2] but at small \B±\ < 0.1 T. All curves are sym- 
metric in respect to the magnetic field direction, despite 
a small offset of -Boffset — —2 mT caused by a residual 
background magnetic field. As expected, the KDC for 
N = 2 (inset) as well as the co-tunneling differential 
conductance for N = 1 (not shown) decrease monoton- 
ically when the magnetic field increases. Surprisingly, 
for N = 1 the KDC shows a non-monotonic behav- 
ior. A pronounced local minimum at B± ~ -Boffset — 
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FIG. 3: Non-monotonic magnetic field dependence of the 
Kondo effect in a one electron DQD. Measurements of the 
KDC as shown in Fig. [2] but for the small magnetic field 
limit. The inset and the main figure plot data for N — 2 
and N = 1, respectively. One data set (open circles) is ver- 
tically shifted for clarity. The dashed curve describes the 
limit B± <C Bk for N = 1 and identical parameters as the 
corresponding model curve in Fig. [2] Vertical lines mark lo- 
cal minima of the N = 1 KDC. Both datasets are taken at 
very similar conditions. Nevertheless, in these as in all our 
measurements we observe slight differences in the oscillation 
amplitude. 



is followed by a quasi-periodic oscillation with minima 
at \B ± -B offsct \ ~ 0,20, 30, 42 mT (vertical lines in 
Fig. [3]). These oscillations quickly decay with increas- 
ing magnetic field and are convolved with the expected 
decrease of the KDC for B <C Bk (dashed line). 

One can consider the following possible explanations 
for the KDC observed in a DQD for TV = 1 to be a 
nonmonotonic function of B±: Namely, (I, II) the leads, 
(III) nuclear spins, (IV) Aharonov-Bohm (AB) like inter- 
ferences, or (V) the alignment of energy levels of the two 
adjacent QDs. All but the last of these possibilities can 
be ruled out. 

(I) Shubnikov-de Haas oscillations in the leads cannot 
depend on the number of electrons. Thus, in contrast 
to our findings they should manifest identically in both 
cases for TV = 1 and 2. 

(II) Spin orbit (SO) interaction in the leads may re- 
sult in a suppression of the KDC in zero magnetic field 
due to spin entanglement between the electrons in the 
lead and the dot [la ]. The SO interaction, however, can 
not explain oscillations of the KDC and a quasi-periodic 
peak structure. Furthermore, the SO interaction should 
equally affect the dot whether charged with one or two 
electrons. Experimentally, there is no evidence of a con- 
ductance minimum for N = 2 (see inset of Fig. [3]) . We 
conclude, that the influence of SO interaction is negligi- 
ble for the effect we observed. 

(III) In GaAs, ~ 10 5 nuclear spins form an internal 



Overhauscr field B nuc ~ 10 mT applied to the electrons 
on each QD. This field fluctuates on the time scale of 
t/V ~ 10 ms [T7[. In our lock- in measurements the data 
are averaged over a much longer time of ~ 300 ms. Hence, 
the fluctuations of B nuc are unlikely to be responsible for 
the observed oscillations. 

(IV) Interference effects in the orbital motion of an 
electron in a double well potential, which determines 
the DQD, could lead to AB-like oscillations in the am- 
plitude of the tunneling between the two wells. In 
terms of the magnetic flux the period of the oscilla- 
tion is the flux quantum $o = h/e [li| . The overlap 
of the wavefunctions centered in the adjacent wells of a 
DQD is proportional to a relative phase shift (V'ilV'2) oc 
cxp(27ri73S'b/$o) of the classically forbidden region Sb ~ 
Ld (we assume a rectangular barrier of width d, lat- 
eral extension L and height V separating the two wells). 
Therefore, one would expect the period of the oscillations 
to be of the order & /S h which is ~ 0.5 T for our DQD. 
This is far in excess of the typical quasi-period (~ 10 mT) 
of the observed oscillations (Fig. [3]). 

As to AB interferences between different tunneling 
paths, the area enclosed by a possible AB contour can 
be estimated to be Sab ~ Ld (or even smaller due to the 
serial configuration of the DQD). 

We also consider the possibility of Fano-like oscilla- 
tions, which would be associated with a leakage current 
under our gates [l9l |. However, experimentally we ex- 
clude leakage currents due to several measurements (not 
shown here). 

(V) We believe that the observed quasi-periodic os- 
cillations can be attributed to the magnetic field effect 
on the alignment of the energy levels in the two adja- 
cent QDs. There is no reason to expect that the DQD 
structure is perfectly symmetric and the single electron 
cigen-encrgies £1.2 in the two QDs are exactly identi- 
cal. The transmission through the barrier correspond- 
ing to a one-dimensional motion of an electron is deter- 
mined by the overlap of the wave functions (V'ilV'2) ~ 
sinh(K_)e~ K +/K_, where k± = (k± ± K%)/2 and k 12 = 
y /, 2md 2 (V — £1,2)/^- The tunneling rate reaches its max- 
imum at resonance for n\ = k%. 

Since the eigen-energies and, consequently, the pa- 
rameters Ki,2 depend on the magnetic field, the latter 
can be used to fine tune the tunneling rate. What 
is the magnetic field that can compensate a mismatch 
Ae = \e\ — e!j I between the ground state eigen-energies of 
the two wells? Using the ld Schrodinger equation for a 
rectangular (or parabolic) double well potential one ob- 
tains Huj c ~ V ' £ i e 2\/ '^£ /W ', where W is the energy dif- 
ference between the two local minima of the double well 
potential and lj c — \eB/m\ is the cyclotron frequency. 

At B± = the energy mismatch Ae is finite. It takes 
about Bq ~ ±12 mT to align the ground states, which 
corresponds to the first KDC maximum. The suppres- 
sion of the mismatch by the magnetic field explains the 
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pronounced minimum at B± = 0. This behavior is ob- 
viously symmetric in respect to the sign of the magnetic 
field. Note, that the characteristic magnetic field B can 
be very small due to the factor ^Ae/W, roughly esti- 
mated to be ~ 10 -3 for our setup. The interdot tunneling 
rate is unaffected by thermal line broadening as long as 
it corresponds to low frequency noise allowing adiabatic 
alignment of the energy levels in both QDs. 

The remaining KDC maxima at slightly larger mag- 
netic fields probably correspond to alignments of excited 
energy states in the two QDs. This implies the impor- 
tance of co-tunneling processes, which are indeed strong 
(compare Fig. [2(c)). Moreover, in order to explain the 
observed quasi-periodic magnetic field dependence, the 
level structures in the two QDs should differ due to some 
anisotropy. From the number of observed KDC maxima 
we conclude that at least three excited states are involved 
in the co-tunneling. 

Our model is consistent with the missing KDC oscil- 
lations for the doubly occupied DQD (N = 2). Indeed, 
for N = 2 and a symmetric charge distribution with one 
electron in each dot, transport is determined by the sin- 
glet and triplet states. The symmetric charge distribu- 
tion allows enhanced elastic co-tunneling decreasing the 
dependence on the misalignment of the single dot ground 
states. 

In conclusion, we here presented measurements of the 
Kondo effect on a DQD charged by only one or two elec- 
trons. We demonstrate control of the resonant tunnel- 
ing in the one electron case by means of a magnetic field 
which appears to be surprisingly small. A non-monotonic 
magnetic field dependence of the KDC is attributed to 
the anisotropy of the DQD. The magnetic field fine tunes 
the alignment of energy levels in the adjacent QDs, mod- 
ifying the interdot tunnel splitting. Hence, the magnetic 
field provides an extremely sensitive tool to detect and 
control the anisotropy of a single electron DQD. 

Our double dot with strong coupling to the leads is not 
a perfect system for a qubit. However, the influence of 
a very small misalignment of the single dot energy levels 
onto the tunnel probability of the electron will remain 
even for a double dot weakly coupled to the leads, which 
can be used e.g. as a charge qubit [2(], 21 1. The ob- 
served variations of the interdot tunnel splitting at very 
small magnetic fields have to be taken into account for 
the design of semiconductor nano-devices in the field of 
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